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Abstract. For a given external loading on a structure we consider the 
optimal stresses. Ignoring the material properties the structure may 
have, we look for the distribution of internal forces or stresses that is in 
equilibrium with the external loading and whose maximal component is 
the least. We present an expression for this optimal value in terms of 
the external loading and the matrix relating the external degrees of free- 
dom and the internal degrees of freedom. The implementation to finite 
element models consisting of elements of uniform stress distributions is 
presented. Finally, we give an example of stress optimization for of a 
two-element model of a cylinder under external traction. 



1. INTRODUCTION 

This paper presents an analysis of optimal stresses in structures under 
given loadings. In P3E1IB], optimal stress distributions for continuous bod- 
ies were considered. Although the problem for a continuous body is more 
difficult mathematically, the corresponding analysis for a structure having a 
finite number of degrees of freedom is more relevant for engineering applica- 
tions. 

Prom the point of view of statics, engineering structures — starting from 
simple trusses all the way to finite element models used for stress analysis of 
continuous bodies, and evidently, continuous models of bodies — are predom- 
inantly statically indeterminate. Mathematically, this means that we have 
more unknown parameters describing the stress distribution in the structure 
under consideration than equilibrium equations. This mathematical problem 
is solved usually by the introduction of constitutive relations and by coupling 
the statics problem with kinematics. 

The work presented here takes a different approach to statically indetermi- 
nate problems. Remaining within the framework of statics, we do not specify 
any constitutive relations and look for the values of the unknown components 
that satisfy the equilibrium conditions and for which the maximal compo- 
nent is the least. Specifically, the problem may be stated as follows. Let 
fd, d = 1,...,D, be the components of the known external loading vector 
on the structure where D is the number of degrees of freedom the structure 
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has, and let cp n , n = 1, . . . , N, be the components of the unknown vector of 
internal forces — stress-like entities. As we consider statically indeterminate 
problems, N > D. The equations of equilibrium will be of the form 

A T &) = f, or A T dn <p n = fd, (1-1) 

where we use the summation convention and the reason we write the matrix 
as A T rather than simply A will be made clear below. Thus, letting c/9 max = 
max n \(p n \ , we are looking for 

S° pt = min{<y9 max } , (1.2) 

where the minimum is taken over all <p satisfying A T (ip) = f. 
Our basic result states that 



■iff* = max (1.3) 



where the maximum is taken over all global virtual displacement vectors 
w = (u>i, . . . , wd) of the structure, so fdWd is the virtual work performed by 
the external force vector. 

A related quantity that we consider is the stress sensitivity of the structure 
defined as follows. Assuming that the internal forces have the same physical 
dimension as the external forces (dimensions of forces or forces divided by 
area), consider the ratio 

gopt 



'V . ... • (1-4) 



max d \f d \ 

Thus, Kf measures the sensitivity of the structure to the external force /. 
Next, we let the external force vary and we look for the worst possible ratio. 
The stress sensitivity of the structure is defined as 

r s° pt } 

K = max Kf = max < — —-. > . (1-5) 

/ / [maxdl/dlj 

It is shown in Section l3~5l that 

ts £—id \ w d\ /1 r s 

K = max — . — A r. (1.6) 

w lun \ A ndW d \ ' ' 

We emphasize that K is a geometric, kinematic property of the structure, 
i.e.., independent of material properties, loading conditions, etc. 

The paper is outlined as follows. We start with the notation and basic facts 
regarding statically indeterminate structures. Then, we prove the results 
stated above. For the internal force vector ip, the value (p ma x is represented 
as a norm, specifically, the dual of the norm ||x|| = \Xn\ that we use 
for internal displacements \- The basic tool we use is the norm preserving 
extension of functionals (the simplified, finite dimensional case of the Hahn- 
Banach theorem). Next, we present some details regarding the application of 
the method to the case of finite element models consisting of elements having 
uniform stress distributions. Finally, in a way of example, we consider the 
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Figure 2.1. A SIMPLE INDETERMINATE STRUCTURE 

case of a two elements model of a thick cylinder under external symmetric 
loading. 

2. STATICALLY INDETERMINATE STRUCTURES 

An elementary example for the type of structures we consider is shown in 
Fig. 1. Our method applies to a lot more complicated structures including 
a large variety of finite element models. 

2.1. Kinematics. The structure is assumed to have D degrees of freedom. 
This means that we have a D- dimensional vector space W containing the 
external infinitesimal virtual displacements (generalized velocities) that are 
compatible with the displacements boundary conditions at the supports and 
the various constraints implied by the structural connections. A generic 
virtual displacement in W will be denoted by w. Thus in our model example 
of Fig. 1, the structure has two degrees of freedom and a generic virtual 
displacement is of the form w = {w\ 1 W2). Each external degree of freedom 
induces a base vector in W. 

Next we consider the space S containing noncompatible, or internal, infin- 
itesimal virtual displacements of the structure. Such an internal deformation 
field of our model example is shown in Fig. 1. Here, the constraints of the 
structural interconnections of the various structural elements are not kept 
(the connection at the bottom joint in the figure). In particular, constant 
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strains within structural elements may be represented as internal virtual dis- 
placement fields. A generic internal virtual displacement field will be denoted 
as X- Clearly, internal deformations have more degrees of freedom than the 
external ones. We assume formally that N, the dimension of 5, is strictly 
larger than D. The various internal degrees of freedom induce base vectors 
in S. (Returning to our model example, we note that the internal degrees 
of freedom, i.e, the base vectors in 5, may be unit changes in lengths of the 
bars or unit axial strains in the bars. It will be convenient also to define 
for a uniform strain structural member, such as a bar in the example, base 
vectors in S consisting of unit strain components in the element multiplied 
by its volume as in Section HJ) 

An important role in the kinematics of the structure is played by the 
interpolation mapping 

A: W — >S. (2.1) 

This mapping associates an internal deformation vector with every external 
deformation. Clearly, as N > D 7 not all internal deformations may be ob- 
tained as images of external deformations under A. We further note that in 
case the structure is not supported distinct displacements fields that differ 
by a rigid displacement field induce the same strain field. While unsup- 
ported bodies may be considered following the methods of j2], we simplify 
the analysis and assume that the supports prevent such rigid displacement 
fields. Thus, we assume mathematically that the interpolation mapping is 
one-to-one, so the matrix of A is of full rank D. 

2.2. Statics. An external force / performs virtual work (power) for various 
external displacements. Denoting by fd the component of the force dual to 
the degree of freedom Wd, the virtual work may be written as 

fW = fdw d - (2.2) 

In other words, we regard an external force as a linear functional 

f:W — >R (2.3) 

and the collection of all external forces is the dual space W* containing all 
real valued linear mappings defined on W. 

In analogy, an internal force ip performs virtual work for virtual internal 
deformation fields. Denoting the component of the internal force <p corre- 
sponding to the component Xn, by tp n , n = 1, . . . , N, we write for the internal 
virtual work performed by ip for the virtual displacement x 

fiX) = fnXn- (2.4) 

Thus, an internal force is a linear mapping 

p:S — ► R, (2.5) 

i.e., cp belongs to the space S* of real valued linear mapping on S. Note 
that in case the component Xn indicates a constant component of the strain 
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in some structural element multiplied by its volume, then ip n indicates the 
corresponding stress component (see Section 0J. 

The principle of virtual work serves as the condition for equilibrium within 
the framework of the structural model. Using the notation introduced above 
it states that 

<p(A(w))=f(w), (2.6) 

for all external vector fields w in W. Using matrix notation where we keep 
the same symbol for a linear mapping (or a vector) and its corresponding 
matrix (or the corresponding column vector) the principle of virtual work is 
written as 

ip T Aw = f T w. (2.7) 

Thus, the principle of virtual work, or equivalently the equilibrium condition, 
may be written in any of the following forms 

<poA = f, A T p = f, A nd ip n = f d , A*(<p) = f, (2.8) 

where in the last equation above we used the dual mapping A* : S* — ► W* 
defined by the condition A*(ip)(w) = (p(A(w)) and whose matrix is the 
transpose of that of A as expected. 

Given an external force /, the equilibrium conditions l|2.8p provide a sys- 
tem of D equations for the N components of the internal force (p. As it 
was assumed that N is strictly larger than D and that A is one-to-one, 
this system of equations cannot determine ip uniquely. In fact, there is an 
(N — D)~ dimensional vector space $ = (A*)^ 1 {/} of solutions to the equi- 
librium problem. 

3. OPTIMAL SOLUTIONS 

For given structures, where the material properties of the various struc- 
tural elements are known, the constitutive relations provide the additional 
information so the internal force vector can be calculated uniquely for any 
given external force /. Here however, we consider the situation where no 
constitutive relations are given a-priori, and among all solutions <p in $, we 
look for the least bound on the maximal component. 

Specifically, for each internal force (p, we set 

IMI°° = V^max = m ax \(f n \ , (3.1) 
n 

and we look for 

S° f pt = min {IMI 00 } = min \ max \ip n \ > , (3.2) 
J tp ip I n ) 

where the minimum is taken over all internal forces (p satisfying A*(tp) = f. 
Thus, for the case were the components of the internal force represent stresses 
in the various structural elements, we are looking for the least bound on the 
discretized approximating stress field. 

In order to evaluate S^ pt directly on the basis of its definition, one would 
have to generate the space of solutions $ and then evaluate the optimal 
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bound in that space. The analysis we present in the sequel will give an 
expression for 5j Pt that does not require the solution of the equilibrium 
equations l|2.8|) . 

3.1. A Solution as an Extension of a Functional. Since the interpola- 
tion mapping A is one-to-one, its inverse 

A' 1 -.Image A — >W (3.3) 

is well defined on its image, a subspace of S. Given an external force /, we 
consider 

f = f o A' 1 : Image A — ► R, (3.4) 

a linear mapping defined on the subspace Image A. Note that the equilibrium 
condition f(w) = ip(A(w)), for a linear functional (p on S may be written as 

l p( x )=foA- 1 ( X ) = f(x) (3.5) 

for all compatible internal fields x in Image A In other words, a solution 
of the equilibrium equations is a linear mapping (p defined on S that agrees 
with / = / o A^ 1 on the subspace Image A Thus, <p is an extension of / to 
the entire space S. 

It is noted that algebraically, given the linear functional / on the subspace, 
it is straightforward to generate an extension of it to the entire space. In 
fact, it is sufficient to show this for the case where N — D = 1, so we have 
to extend / to a space having one more dimension. In the general case 
where N — D is any other finite number, the procedure can be carried out 
inductively adding one dimension at a time. To generate such an extension 
if, one can choose an internal displacement vector xi that does not belong 
to Imaged and give an arbitrary value to c = f(xi)- Any vector \ m the 
larger space may be written as a linear combination 

X = Xo + axi (3-6) 

for a compatible internal virtual displacement xo m Imaged and a real 
number a. Thus, for any extension tp we have 

V(x) = ^(Xo + axi) 

= <p(Xo) + apixi) 

= /(xo) + ac. (3.7) 

3.2. Norms. Recalling that we are looking for a solution of the equilibrium 
equations that minimizes {{(fW 00 = max n \cp n \, we mention a number of useful 
properties of this norm on the space of internal forces. Consider the norm 
on the space S given by 

llxl|i = £)lX"l- (3-8) 

n 
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Then, the following holds 

II ll°0 \fi.X)\ /q Ci\ 

99 = max -r — n — . (3.9) 

* llxlli ' 

This relation between the norm for internal forces and the norm for internal 
virtual displacements is all we need. In fact, the following analysis applies 
to other criteria for optimization of the internal forces, i.e., criteria given 
by other norms say \\<p\\. To do this, one should determine the norm on 
the space of internal displacements to which \\ip\\ is dual, i.e., determine the 
norm ||x|| such that 

II II |V(X)1 

w =max— — — . (3.10) 

* llxll 

Thus, in the sequel we use a generic norm ||x|| for internal displacements 
and the corresponding dual norm \\ip\\ for internal forces satisfying Equation 
(j3.1()JI . In fact, once the norm \\x\\ is established as the one for which our 
optimality condition is a dual norm, the norm on the space of internal dis- 
placement fields (rather than the one on the space of internal forces) plays 
the central role as will be seen below. It is noted that if the optimality cri- 
terion is given in terms of a norm \\ip\\, the associated norm on the space of 
internal displacements can be found by (see gj p. 186]) 

||X||= max M#. (3.11) 
<p \\(p\\ 

For external virtual displacements we may also consider the norm 



d 



and use the dual norm 



CO 



max- — - — = max (3.13) 

u> \\w\\i d 



for external force vectors. 

Contrary to the previous paragraph, we will find it useful in the analysis 
of optimal internal forces below to use for external virtual displacements the 
norm 

HI = PMU, (3-14) 

where on the right we use the norm on S. The corresponding dual norm for 
external forces is therefore 



(3.15) 
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3.3. Optimal Extensions. Considering the linear functional / = foA 1 : Imag 
M, we may evaluate its dual norm relative to the norm l|3.8p on Image A 
Thus, 



max 

XGlmage A 



fix) 



rali 



(3.16) 



where it is noted that the maximum is evaluated for all x m Imaged (and 
not the entire space S). The Hahn-Banach theorem of functional analysis 
states that there is a linear functional (/?j Pt : S — > M (i.e., defined on the 
space <S) such that 

^7\x) = fix) 

for all x m Imaged, and 



(3.17) 



opt I 



max 



l/ll 



max 



fix) 



XelmageA X 1 



(3.18) 



In other words, ip extends / without increasing its norm. It was mentioned 
earlier that extending the functional one dimension at a time is simple (see 
Equation l[3.7[) . Appendix A presents the construction for the addition of one 
dimension to the domain without increasing the norm of the linear functional. 
In fact, the Hahn-Banach theorem asserts that this can be done for infinite 
dimensional spaces also and is used in the continuum counterpart of this 
analysis in [H |2| IHj ■ 

We end this subsection by noting that 



max 

XGlrnage A 



fix) 



rail 



l/HI 

whw \\A(w)\\ 



max 



(3.19) 



3.4. The Equation for the Optimum. Returning to the expression for 
the optimum, we note that in general for any internal force <p that extends 
1 



max 



fix) 



XSlmageA ||x|| 

\vix)\ 



max 



XGlmageA ||x| 

\<pix)\ 



< max 



xes Hxl 

Opt II 



(3.20) 

(by the principle of virtual work) (3-21) 
<P\\°°- (3-22) 



Thus, since for - op1 n i |x 
we have 



\f\\°° = \\f\\, as in Eqs. (l3~T8l and llTTfll . 



iv/ pt ir 



mm \\(p\\ 



(3.23) 
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where the minimum is taken over all extensions <p of /. Since an extension 
if of / satisfies the equilibrium condition as in Eq. (|3.5|l , we conclude that 



mm < max 

f=A*( V ) L n 



.opt 1 1 CO 



It follows that 

max — — — . (3.25) 



•.opt \f(w)\ \fdW d \ 
•> / = max ,, —. ... = max — 

' tuGW ||A(«>)|| rtW J2n \ A ndW d \ 



3.5. The Stress Sensitivity. We now derive the expression for the stress 
sensitivity of the structure. Recalling its definition in Eq. (|1.5|) 



K = max 



<^opt 



/ l^maxdl/dlj 

we substitute the expression for the optimum internal force to get 

^ = max{ \- (max ^ 1 ^ 1 , U (3.26) 

™ eVV I l^nd^dl I / max d \f d \ J J 
We may use now Eq. (|3 .11^ and arrive at 

K = max rP?} Wd \ ; • (3.28) 

™ eW En lAid^d! 



4. APPLICATION TO FINITE ELEMENTS 

A typical situation where one would like to apply the foregoing analysis is 
a finite element model of a continuous body. Thus, we briefly describe here 
some additional details for the simple situation of a finite element model 
where it is assumed that the stress is uniform within each element. We 
do not consider here the question of approximation of the solution to the 
continuum problem by finite elements and take the finite element model as 
a given structure. 

Let L be the number of elements, e; the l-th element, and a\j the com- 
ponents of the uniform stress in that element. We want the collection of <r| - 
for the various elements I and various components i,j to be the components 
of our internal force vector. Thus, the index n is replaced by the collection 
of 3 indices I, cp^ = crv, and 

(4.1) 

is the quantity that we want to minimize. 

The internal degrees of freedom v;-.- should be chosen such that Y\ , • ctL yL- 
is the virtual work of the internal force for the given internal displacements. 



\ip\\ = max 

i.n.l 
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Writing 



U = J aijEij dV = Y, o\j J *j dV, 



(4.2) 



where denotes the linear strain field, one realizes that the internal degrees 
of freedom are given by 



&, = j 



Eij dV. 



(4.3) 



For uniform strain elements x\j = ^l £ lj ( n0 sum 011 0> wnere £ \j are the 
components of the uniform strain in the element e; and Vi is its volume. The 
norm ||x|| is therefore 



iixii = E 



£ij dV 



and for uniform strain elements 



iixii = J2\M Vl - 



(4.4) 



(4.5) 



It is a standard procedure in the construction of finite element models to 
use an array E\- d such that 



Xij 



Eij dV = E ijd w d , 



(4.6) 



where the array E\- d is generated using the relation between the external 
degrees of freedom and the displacements at the nodes of the the element 
ei, using the shape functions to obtain the displacement field within the 
element, and using the discrete forms of the differentiation and integration 
operators in order to obtain the strain components and their integrals. Thus, 
the array E\- d replaces the matrix A n d in our expression for the optimum. 
We conclude that 



r<opt 



I opt I 
\<Pf I 



mm 



max 



a. 



max 

wew 



\fdW d \ 



Emj E \jd w d 



(4.7) 



5. EXAMPLE 

We consider the following 2-dimensional example. The structure is a finite 
element model of a cylinder of inner radius r\ and outer radius r% under 
external normal traction p at the outer boundary (see Fig. 2). 

The finite element model consists of two uniform stress elements e\ and 
e-2, corresponding to the regions r\ ^ r ^ r2, and ^ r ^ r%, respectively, 
where T2 = (ri + r^)/2. Due to the cylindrical symmetry the problem has 
three degrees of freedom and a typical external virtual displacement is of 
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the form w = (wi,W2,Ws), where Wd is the radial displacement at r^. The 
corresponding components of an external force / are given as 

fd = 2-KrdPd-, (no summation), d=l, 2, 3, (5.1) 

where, pd is the applied load at r^. Thus, for the case under consideration 
/ = (0, 0, 2nr3p). The space of internal virtual displacements S will be 4- 
dimensional and an internal displacement will be of the form xi = 5*i(e r )i, 
X2 = Si (eg) i, X3 = S 2 {£ r h, Xi = S2{ee)2-, where 5/ is the area of the Z-th 
element and (e r )i, (ee)i are the uniform strain components in e/. The values 
of the strain components within the elements are approximated as 

_ wi +l - Wl _ Wl + w i+1 

n+i - n 2m i 

where mi = (r; + n + i)/2 denotes the mean radius of the Z-th element. The 
matrix A is easily calculated to give 
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Instead of the optimal stress, it will be convenient to determine 

1 • En \ A ndW d \ . . 

— mm — j- : . (5.4) 



S°pt w \f d w d \ 

Since both numerator and denominator are homogeneous in the vector w, we 
normalize it by conducting the search for the minimum over all displacement 
vectors satisfying 

fdW d = fow 3 = 1 (5.5) 

so 

W3 = ' (5 ' 6) 

and we have to minimize _ 

\A nd w d \ (5.7) 

n 

over all values of (wi,W2) subject to the condition 1)5. 6j) . Writing the sum 
above explicitly and using the fact that A 3 i = A 32 = 0, we have 

— — = min{P+ \A 3d w d \ + \A Ad w d \} , (5.8) 

where 

P = min {\AnWi + Ai 2 w 2 \ + \A 2 \Wi + ^22^2!} • (5.9) 

As a function of w\, the expression in the curly brackets above is piecewise 
affine and attains its minimum at some point where two adjacent line seg- 
ments of its graph meet — at some value of w\ where one of the absolute value 
terms vanishes. This gives 

p _ |u> 2 | |^4ll^22 ~ ^12^21 1 _ \w 2 \ \ A U A 2 2 - ^12^21 1 .„x 

max{|Aii|,|A 2 i|} \A U \ ' ' 

as evidently \Au \ > \A 2 i\. Substituting the values of the various components 
of the matrix we obtain 

P = 2vr(r 2 - n) \w 2 \ (5.11) 

We can turn back to the minimization of Eq. 1)5.8)1 where now the minimum 
is attained at a value of w 2 where one of the two absolute value terms vanishes 
or at the value where P vanishes (so w 2 = 0). Setting 

Q = P+ \A 3d w d \ + \A 4d w d \ , (5.12) 

for the case A 32 w 2 + A 33 w 3 = 0, we have 

Q = Qi = j?i ttA □ n + D 4 , (5.13) 

\A 32 \ I max { I An [ , 1-421 1 } J 
where D x = A n A 22 - Ai 2 A 2 i, and D 2 = A 33 A 42 - ^43^32 = 2ttS 2 are the 
determinants of the two submatrices. It follows that 

Q 1 = 2 7 r(r 3 -r 1 )K|. (5.14) 
For the case w 2 = 0, we have 

Q = Q 2 = {\A 33 \ + \A 43 \) \w 3 \ = 2vrr 3 \w 3 \ , (5.15) 
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and for the case A i2 w 2 + ^43^3 = 0, we have 




^43 1 D 1 D 2 



A 42 \max{\A 11 \ ,\A 2 i\} \A i2 



= 27r(r 3 + 2r 2 - n) \w s 



(5.16) 



It follows that 



1 



min{Qi, Q 2 , Q3} = Qi 



(5.17) 



^opt 



and substituting the normalization condition on 103, we finally obtain 



It is noted that this value corresponds to a uniform value of uq that will 
balance the external loading for one half of the cylinder. 

This result for a common engineering problem could have been possibly 
obtained by direct analysis of the hollow cylinder. However, this example 
was aimed at demonstrating the strength of the proposed procedure and 
the ability to implement it in real life structures. Clearly, aside from the 
practical significance of knowing the optimal stress that may develop in a 
structure under given loading conditions, this procedure can be used to assess 
the optimality of a given standard design. Thus, comparing of the maximal 
stresses developing in a proposed design with the optimal value obtained by 
application of the formulation outlined above, one can estimate how much 
the design can be improved. 
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Appendix A. Appendix: Norm Preserving Extensions 
Let Wo C W be a vector subspace and Fq : Wq — ► K a linear functional 



P- 



r 3 - n 
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for all wq € Wq. 



(A.l) 
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Let W\ be an element of W — Wo and let W\ be the vector space spanned 
by Wo and w\, i.e., 

Wi = {w + awi; a G R, w G W } . 
We want to extend Fo to a functional 

Fi : Wi — ► R, 

such that 

F 1 (w )=F (w ), for all w eW (A.2) 

Fi(«j) < HFolllHI, for all w G Wi- (A. 3) 
Assuming Fi satisfies the above requirements, then, 

F x (u;o+^i) = Fi(w ) + Fi(wi) (A.4) 

ll-^ollll^o + will > F (^o) + Fl(wi) using II Fi|| = ||F II . (A. 5) 

Thus, 

Fi(toi) ^ ||F ||||w;o +toi|| - F (w ) for all wq G Wo- (A. 6) 
Similarly, 

Fi(-wo-wi) = -i*i(w )--Pi(u>i) (A.7) 

||F |||| - too -will > -Fo(too) -Fi(tyi) (A.8) 

||F |j|K + HI > -F (ty ) -Fi(toi). (A.9) 

Thus, 

FiOi) ^ -||Fo||||ioq + toi|| -F (w ) for all w G Wo- (A. 10) 

We now show that the necessary conditions of Eqs. (|A.6|) and (|A. 10|1 are 
also sufficient. For any w = wo + aw\ € Wi we have 

Fi(u>o + awi) = aF (w /a) + aFi(u>i), (A. 11) 

using F 1 (w ) = F (wo). 

Now if a > 0, we may use Eq. I|A,6|1 for wo/a in the second term on the 
right of Eq. 1|A,11|) to get 

Fi(w + awi) ^ aF (wo/a) + a (||F || \\w /a + w\ \\ - F (w /a)) 
= ||Fo||||ioo + awi\\. 
Alternatively, for a < 0, Equation (|A.10|) may be rewritten as 

aFi(u>i) ^ — a||Fo \\\\wo + w;i|| — aFo(wo) for all wo G Wo- 
When we substitute this for wo/a in Equation 1|A,11|) we obtain 

Fi(u>o + awi) < aF (w /a) - a||F || \\w /a + wi\\ - aF (w /a), 
= -a\\F \\\\wo/a + ti>i||, 
= ||Fo||||wo + awi\\, 
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where we used the fact that (—a) > in the last line above. This completes 
the proof. 
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